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1  INTRODUCTION 

Particle  beam  lithography  depends  upon  the  complex  interaction  of  the 
incident  particle  with  the  polymer  resist  used.  Irradiation  of  such  a  polymer 
by  particles  causes  either  polymer  chain  scission  (positive  resist)  or  polymer 
chain  cross-linking  (negative  resist)  giving  rise  to  areas  of  different 
solubilities.  The  change  in  this  property  can  be  used  in  a  differential 
solubility  process  using  an  appropriate  organic  solvent  to  remove  more  soluble 
areas  leaving  less  soluble  areas  largely  unaffected.  This  creates  a  pattern 
on  the  substrate  that  can  be  used  to  mask  subsequent  processing  steps. 
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A  detailed  understanding  of  the  particle-target  interactions  is  fundamental 
to  being  able  to  produce  very  small  features  routinely.  In  particular  it  is  very 
important  to  be  able  to  understand  and  to  compensate  for  the  inevitable  proximity 
effect  that  is  unique  to  electron  beam  lithography.  This  proximity  phenomenon 
can  be  characterised  as  two  similar,  yet  distinct,  effects: 

i.  Intra-proximity  effect,  which  is  caused  by  forward  scattering  of  the 
incident  electrons  resulting  in  a  reduction  of  exposure  dose  (charge/ 
area)  and  an  effective  broadening  of  the  incident  beam. 

ii.  Inter-proximity  effect,  caused  by  backscattered  electrons  re-emerging 
from  the  substrate  and  interfering  with  the  exposure  of  that  feature 
and  of  neighbouring  features. 

Ion  beam  lithography  is  attractive  because  the  much  heavier  ions  give  up  their 
energy  very  rapidly  and  undergo  very  little  scattering.  Hence,  the  resist  has 
a  greater  effective  sensitivity  to  exposure  by  ions  and  there  are  no  proximity 
problems  to  limit  the  resolution  of  the  process.  However,  it  is  a  relatively 
new  technology  with  which  little  has  yet  been  achieved  and  therefore  much  less 
information  is  available  about  the  process  and  its  performance. 

To  help  study  and  understand  the  exposure  processes  of  particle  beam 
lithography  a<qomputer  modelling  system  has  been  developed.  It  was  decided  to 
develop  a  new  suite  of  programmes  in  order  to  provide  a  basis  for  the  investi¬ 
gation  of  the  assumptions  and  approximations  made  in  the  various  models  and  to 
enable  the  rapid  modification  of  the  system  to  tackle  different  cases  as  these 
arise  in  the  complementary  experimental  programme.  At  the  heart  of  the  system 
is  a  Monte  Carlo  simulation  programme.  This  uses  random  numbers  to  return 
values  for  the  scattering  angles  and  step  length  between  interactions  with  the 
target  nuclei  in  order  to  simulate  the  particle  trajectories.  The  programme 
stores  the  energy  lost  by  those  particles  as  they  travel  through  the  resist 
layer  to  build  up  an  exposure  profile  for  the  process. 

The  study  of  electron  scattering  within  a  material  target  using  a  Monte 
Carlo  technique  has  been  investigated  in  detail  by  Bishop'1 ',  using  the 
screened  Rutherford  interaction  cross  section,  Bethe's  continuous-slowing-down 
approximation  and  Lewis'  multiple  scattering  theory  to  calculate  X-ray  produc¬ 
tion  by  an  electron  probe.  Adopting  the  same  approach,  Shimizu  and  Muratat?) 
reported  studies  on  the  resolving  power  and  contrast  of  an  SEM  concentrating 
on  the  backscattered  electrons.  A  change  to  using  the  more  powerful  single 
scattering  model  was  made  by  Murata  et  alO)  and  followed  by  Kyser  and  Murata^' 
in  their  study  of  microprobe  analysis  of  thin  films  on  substrates.  The  first 
application  of  this  technique  to  electron  beam  lithography  was  made  by  Shimizu 
et  al'^'  and  Kyser  and  Murata^)  which  included  studies  of  energy  deposition 
into  the  resist  layer. 

Further  to  these  studies  Adesida  et  al^  have  reported  results  using  an 
extension  of  Gryzinski's  excitation  function  as  a  replacement  for  Bethe's 
continuous-slowing-down  approximation  in  an  attempt  to  simulate  the  exposure 
mechanisms  more  directly.  Also  the  effect  of  secondary  electrons  within  the 
exposure  process  has  been  investigated  by  Mjrata  et  al(®),  but  these  are 
beyond  the  scope  of  the  work  presented  here. 
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Computer  studies  of  ion  scattering  within  solid  targets  have  been  reported 
by  many  authors.  Included  amongst  these  are  several  dedicated  Monte  Carlo 
programmes  such  as  MARLOWE  by  Robinson  and  Torrens^).  TRIM  by  Biersack  and 
HaggmarkOO)  and  PIBER  by  Adesida  and  Karapiperis^  .  MARLOWE  was  initially 
developed  to  study  atomic-displacement  cascades  with  a  binary-collision  approxi¬ 
mation,  assuming  nuclear  scattering  governed  by  the  Moliere  potential.  TRIM 
also  uses  Moliere's  potential  but  incorporates  a  more  efficient  analytical 
method  for  determining  scattering  angles  to  study  ion  transport  in  amorphous 
targets.  PIBER  was  developed  to  study  ion  exposure  of  resists  for  lithographic 
purposes  and  uses  the  experimentally  derived  universal  scattering  cross-section 
of  Kalbitzer  and  0etzmann(13)  and  like  the  previous  two  uses  the  electronic 
stopping  formulae  of  Lindhard  et  al^)  ancj  gethe. 

In  the  programmes  presented  here  the  scattering  model  of  Kyser  and  Murata^&^ 
is  used,  with  alterations,  to  simulate  the  scattering  of  electrons  and  Bethe's 
continuous-slowing-dovn  approximation  is  assumed  to  supply  a  rate  of  loss  of 
energy  by  the  electrons.  A  method  using  Kalbitzer  and  Oetzmann's  universal 
scattering  function,  akin  to  that  of  Adesida  and  Karapiperis^ '*12) f  is 
employed  for  the  ion  case  and  electronic  stopping  is  assumed  by  Bethe's  formulae 
for  high  energy  ions  and  by  the  LSS^)  theory  in  the  low  energy  regime. 

The  programmes  have  been  designed  to  be  flexible  to  allow  the  user  to  define 
an  incident  beam  voltage  incident  upon  an  amorphous  resist  structure,  of  up  to 
three  layers,  on  top  of  an  amorphous  substrate  of  semi-infinite  thickness.  They 
were  developed  to  study  exposure  processes  rather  than  to  be  a  detailed  examina¬ 
tion  of  particle  scattering.  It  is  intended  that  the  simulation  progranmes 
described  within  will  be  expanded  to  cover  the  whole  lithographic  process  providing 
a  comprehensive  simulation  package. 

2  THEORY 

The  scattering  model  of  Kyser  and  Murata^*^  is  used  in  the  computer 
simulation  of  electron  beam  exposure.  For  the  ion  beam  case  the  approach  of 
Anderson  and  Ziegler^'  is  followed,  ie  the  theory  of  Lindhard  et  al^)  £s 
assumed  for  low  ion  energy  levels  and  at  high  energy  levels  the  Bohr-Bethe 
equation  is  applied,  using  the  interpolation  scheme  of  Verelas  and  Biersack 
for  the  intermediate  energy  range. 

2.1  ELECTRON  SCATTERING  MODEL 

It  is  assumed  that  electrons  suffer  Rutherford  scattering  within  a 
binary  collision  process  with  a  target  atomic  nuclei.  The  Rutherford 
differential  scattering  cross  section  can  be  expressed  as 
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where  Z^  is  the  atomic  number  of  the  scattering  nucleus,  e  is  the  electron 
charge,  c0  is  the  permittivity  of  free  space,  E  is  the  energy  of  the 
electron,  6  is  the  scattering  angle  and  a£  is  the  screening  parameter. 
Strictly  equation  (1)  is  not  Rutherford's  original  differential  scattering 
cross  section,  the  Zi(Zj+1)  factor  replaces  the  usual  Z?  factor  in  the 


numerator  to  allow  for  inelastic  electron-electron  scattering  and  a  screen¬ 
ing  parameter,  oj,  is  introduced  to  account  for  the  electronic  screening  of 
the  target  nucleus  by  its  orbital  electrons,  o^  is  given  by  Cosslet  and 
Thomas as 
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where  Xe  is  the  wavelength  of  the  electron  and  aQ  is  the  Bohr  radius  of  the 
Hydrogen  atom. 


The  total  scattering  cross  section,  o£,  is  achieved  by  integrating 
equation  (t)  between  the  limits  0-0  to  0-n,  ie 
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2.2  ENERGY  LOSS  BY  ELECTRONS 


Bethe's  continuous  slowing  down  approximation  is  assumed  to  provide  a 
rate  of  loss  of  energy  by  the  incident  electron  as  it  travels  through  the 
target  material.  The  equation  for  this  is 
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where  y  -  1.166,  Na  is  Avogadro's  number,  p  is  the  density  of  the  target 
material,  is  the  weight  fraction  of  the  i^  constituent  atomic  species 
and  is  its  atomic  mass.  Ji  is  the  mean  ionisation  energy  of  the  nuclei. 
The  values  for  Ji  are  given,  for  j£ < 13,  by 
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Otherwise  by 
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2.3  ION  SCATTERING  MODEL 

Ions  lose  their  energy  by  two  major  processes  due,  on  the  one  hand,  to 
elastic  Coulomb  interactions  between  the  nuclear  charges  of  the  ion  and 
target  nucleus  and,  on  the  other,  to  inelastic  interactions  of  the  incident 
ion  with  the  electrons  within  the  target.  At  low  energies  nuclear  stopping 
dominates  the  interaction  and  is  responsible  for  most  of  the  angular 
dispersion  of  the  incident  ions.  At  higher  energies  electronic  collisions 
are  store  important.  The  theory  behind  these  processes  is  much  less  well 
formulated  than  in  the  electron  case,  especially  within  the  intermediate 
energy  range  in  which  we  are  working. 


2.4  ELECTRONIC  STOPPING 

In  the  high  energy  regime  but  at  non-relativist ic  ion  velocities,  the 
electronic  stopping  power  of  the  ion  is  adequately  described  by  the  Bohr- 
Bethe  equation 
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where  Zj  is  the  atomic  number  of  the  incident  ion,  N  is  the  number  density  of 
target  nuclei.  Me  is  the  mass  of  an  electron  and  B  is  Bethe's  stopping  number 

/2M ey2\ 

B  -  Z2  lnl— j — 1  .  (7) 

Here  1  is  the  average  excitation  energy  of  the  target  electrons  and  is  given 
in  equations  (5a)  and  (5b) .  Z2  is  the  atomic  number  of  the  scattering 
nuclei. 

In  the  low  energy  regime,  for  velocities  up  to  e2/fi  the  LSS  theory 

can  be  applied  to  obtain  a  rate  of  loss  of  energy  by  the  incident  ions 
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where  VQ  is  the  Bohr  velocity  and  €e  is  a  dimensionless  constant  of  the 
order  of  zj/6. 

It  is  convenient  to  define  a  reduced  energy 


e  -  EaM2(4irco)/(Z1Z2e2(M1  +  M2)) 
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in  which  Mj  and  M2  are  the  respective  masses  of  the  incident  ion  and  target 
nucleus  and  a  is  a  Thomas-Fermi  constant  given  by 

a  -  0.885a  (Z2/3  +  Z2/3)_i  (10) 

01  i 

and  a  reduced  range 

Pr  -  RNM2a24rM1/(M1 +M2)2  (11) 


where  R  is  the  real  range.  Using  these  parameters  the  LSS  theory  expresses 
a  rate  of  loss  of  energy  due  to  electronic  stopping  as 
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Following  Anderson  and  Zeigler^1^  the  interpolation  scheme  of 
Verelas  and  Biersack 


(14) 


is  used  to  bridge  the  gap  between  the  high  and  low  energy  regimes. 


2.5  NUCLEAR  STOPPING  AND  SCATTERING  CROSS  SECTION 


The  LSS  theory  assumes  the  Thomas-Fermi  inter-atomic  potential  to 
provide  a  differential  nuclear  scattering  cross  section  which  can  be 
expressed  as 

do  *  (t)  -  2ira2t_3/2  f(t*)  dt  (15) 


where 


and  e  is  the  reduced  energy  of  equation  (9).  The  universal  scattering 
function,  f(ti),  can  be  approximated  to 

f(t*)  -  l  t*“"  [1  ♦  (21  t1“m)q]"1/q  .  (17) 

The  values  of  the  constants  £,  m,  q,  are  given  by  Kalibitzier  and  Oetzmann1 
as 

l  «  2.54  m  -  0.25  q  -  0.475 

The  total  scattering  cross  section  is  assumed  to  be  equal  to  N  (after 
Adesida  and  Karapiperis^ ^) ,  ie 


The  nuclear  energy  loss,  AE„,  is  simply  the  energy  transferred  from 
the  incident  ion  to  the  target  nucleus  on  collision 

AEn  -  AM1M2(M1 +M2)"2  E  sin2^  (20) 

where  6  is  the  scattering  angle  in  the  centre  of  mass  system. 
IMPLEMENTATION 


3.1  THE  MONTE  CARLO  METHOD  FOR  CALCULATING  THE  PARTICLE  SCATTERING 
PARAMETERS 


The  probability,  P(6)  dft,  of  an  incident  particle  scattering  into  the 
solid  angle  dft  *  sin(0)  d6dft  is  given  by 

P(6)dfi  -  (o'(e)/o)  dO  (21) 


where  a  is  the  total  scattering  cross  section  and  o' (6)  is  the  fractional 
scattering  cross  section,  hence,  for  the  electron  case 


(22) 


where  Rj  is  a  uniformly  distributed  random  number  between  0  and  1 .  The  ion 
beam  case  is  less  straightforward  to  implement.  The  integral  of  equation 
(15)  cannot  be  directly  evaluated  and  therefore  has  to  be  calculated 
numerically.  A  look-up  table  of  I(t)  versus  t  is  generated  from  the 
integral 
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at  the  start  of  the  computer  programme  which  can  then  be  searched  to  locate 
a  value  for  t  from  which  the  scattering  angle,  6,  can  be  deduced: 
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The  azimuthal  angle  is  simply  given  by 
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where  R2  is  a  uniformly  distributed  random  number  between  0  and  1 . 

The  step  length  between  collisions,  S£,  is  calculated  using  Poisson 
statistics  which  gives 
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Again  R3  is  a  uniformly  distributed  random  number  between  0  and  1  and  X  is 
the  mean  free  path  of  the  particle.  For  a  mixed  element  target  of  i  elements 
each  with  weight  fraction  C£  and  atomic  mass 


(27) 


where  p  is  the  mass  density  of  the  target.  The  probability,  P^,  that  the 
incident  particle  will  be  scattered  by  atom  i  is  simply  its  fraction  of  the 
sum  of  the  cross  sections,  ie 


and  can  be  decided  upon  by  generating  a  random  number,  R4.  Hence,  if  is 
in  the  range  (0-P^)  the  particle  is  assumed  to  be  scattered  by  an  atom  of 
species  i,  otherwise  if  R  is  in  the  range  (P£  -  1)  scattering  is  assumed  by 
a  different  species.  The  actual  species  is  determined  by  comparing  R4  with 
the  ranges  0  to  Pp  Pj  to  (P1+P2)  etc.  In  this  manner  the  scattering  species 
is  decided  upon  in  a  random  yet  weighted  fashion. 

3.2  SEQUENCE  OF  CALCULATIONS 

A  particle  of  energy  E0  is  incident  at  90  deg  to  the  surface  of  a  semi¬ 
infinite  target.  The  x-y  plane  being  defined  as  that  of  the  surface  with 
the  particle  approaching  along  the  z  axis  travelling  in  a  +ve  z  direction 
(see  Figure  1).  The  step  length  SQ  is  calculated  by  the  Monte  Carlo  method 
described  above,  providing  the  penetration  depth  before  the  first  scattering 
event  at  point  1 .  The  energy  lost  AE0  is  deducted  from  E0  to  give  a  new 
energy  Ej.  The  scattering  angles  8j,  $1  and  step  length  S|  are  then  calcu¬ 
lated  to  arrive  at  point  2.  Again  the  energy  lost  is  deducted  from  Ej . 

This  process  is  repeated  until  the  energy  of  the  incident  particle 
diminishes  to  an  insignificant  level,  or  in  the  case  of  electrons,  that 
their  projected  range  is  too  short  and  they  will  not,  therefore,  re-emerge 
into  the  resist  film. 

The  computer  programme  stores  the  energy  deposited  into  the  resist 
layer  in  a  50  by  50  array  which  corresponds  to  an  imaginary  grid  of  parallel¬ 
epipeds  within  the  resist  film  (see  Figure  2).  The  cross  sectional  area 
of  each  parallelepiped  is  dependent  upon  the  overall  dimensions  of  the 
simulation  but  is  typically  100  Angstroms  square  for  electrons  and  10  Ang¬ 
stroms  square  for  ions.  The  programne  calculates  which  parallelepipeds  the 
particle  passes  through  between  collisions  and  its  path  length  within  the 
respective  elements  to  deduce  the  amount  of  energy  deposited  into  each. 

This  energy  is  then  added  to  the  running  total  stored  in  the  corresponding 
array  element. 

3.3  CONVOLUTION  OF  THE  ENERGY  DEPOSITION  PROFILE 

The  computer  programme  described  above  assumes  that  all  the  particles 
are  incident  at  the  same  point  on  the  surface  of  the  resist.  Obviously, 
it  is  impossible  to  produce  a  system  with  such  a  beam  shape,  hence,  the 
energy  deposition  profile  is  convoluted  with  the  Gaussian  beam  profile  of  a 
real  system  to  simulate  an  energy  deposition  profile  for  a  single  pass 
exposure . 
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Let  E(y,z)  be  the  energy  dissipation  per  unit  volume  at  a  lateral 
distance  y  and  depth  z  from  the  origin  for  the  point  source  and  the 
current  distribution  of  the  Gaussian  beam  be  described  by 
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where  rQ  is  the  radius  of  the  beam.  The  beam  diameter  is  defined  by  the 
10Z-90Z  of  beam  current  width.  Then  the  energy  distribution  for  the 
Gaussian  beam  system  is  given  by 
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RESULTS  AND  DISCUSSION 
4.1  TRAJECTORY  PLOTS 

The  figures  3  to  8  show  the  simulated  trajectories  for  a  point  source 
of  electrons  with  incident  energies  of  10,  20  and  50  keV  in  two  different 
thicknesses  of  PMMA  resist  upon  a  silicon  substrate.  They  clearly  show 
that  the  resolution  of  the  electron  beam  is  severely  limited  by  the  scattering 
of  the  electrons  and  exhibit  the  characteristic  "pear"  shape  of  the  scattering 
The  broadening  of  the  beam  as  it  travels  through  the  material  is  caused  by 
forward  scattering  of  the  electrons  and  is  known  as  the  intra-proximity 
effect.  It  is  noticed  that  this  effect  can  be  greatly  reduced  within  the 
resist  layer  by  increasing  the  incident  beam  energy.  However,  it  is  also 
clear  that  with  an  increase  in  beam  energy  there  comes  an  increase  in 
electron  range  and  consequently  an  increase  in  the  area  affected  by  the  back- 
scattered  electrons.  These  backscattered  electrons  can  interfere  with 
neighbouring  exposures  causing  the  inter-proximity  effect.  These  figures 
also  demonstrate  the  fact  that  resolution  can  be  improved  through  using  a 
thinner  resist  layer,  especially  at  the  lower  beam  energies,  because  the 
resist-substrate  interface  is  brought  up  higher  into  the  neck  of  the 
scattering  "pear"  where  the  intra-proximity  is  less.  This,  then,  can  be  a 
very  useful  method  for  improving  the  resolution  of  electron  beam  machines 
of  limited  beam  voltage. 

The  flexibility  of  the  computer  programme  is  exhibited  by  figures  9 
and  10.  Figure  9  shows  the  simulated  trajectories  of  100  20  keV  electrons 

incident  on  to  0.5  microns  of  PMMA  on  a  Gallium-Arsenide  substrate.  The 
problem  of  increased  scatter  due  to  Gallium-Arsenide  having  a  greater 
density  than  Silicon  is  evident.  PMMA  resist  is  used  here  but  the  programme 
is  capable  of  simulating  any  user  defined  film  on  any  substrate.  Figure  10 
shows  the  simulated  trajectories  of  100  20  keV  electrons  incident  on  to  a 

tri-level  resist  structure  consisting  of  0.5  microns  of  PMMA  on  0.1  microns 
of  Aluminium  on  1.0  microns  of  PMMA  on  a  silicon  substrate.  It  is  noticed 
that  there  is  marginally  less  backscatter  than  in  the  equivalent  one  layer 
case  (Figure  6)  due  to  most  of  the  scattering  occurring  further  away  from 
the  top  of  the  resist  layer. 


Figures  11  to  15  show  the  ion  trajectories  simulated  with  100  keV 
Gallium  into  0.5  microns  of  PMMA  on  Silicon,  200  keV  Silicon  into  0.5 
microns  of  PMMA  on  Silicon,  200  keV  Silicon  into  0.2  microns  of  PMMA  on 
Silicon,  150  keV  Silicon  into  0.5  microns  of  PMMA  on  Silicon,  100  keV 
Silicon  into  0.5  microns  of  PMMA  on  silicon  respectively.  These  two 
species  have  been  simulated  as  they  are  available  from  the  liquid  metal 
sources  to  be  supplied  with  the  RSRE/VG  100  kV  ion  beam  lithography  machine 
(IBL100).  Comparing  these  figures  with  those  for  the  electron  beam  case 
the  potential  for  greater  resolution  is  obvious.  However,  their  extremely 
short  range  may  cause  some  processing  problems.  (This  will  be  covered 
further  in  a memorandum  yet  to  be  published  on  focussed  ion  beam  systems  by 
V  Mifsud.) 

4.2  ENERGY  DEPOSITION  PROFILES 

The  energy  deposition  results  from  the  Monte  Carlo  simulations  of 
electron  beams  have  been  convoluted  with  a  Gaussian  beam  current  distribu¬ 
tion  with  a  10%  to  90%  width  of  0.2  microns  to  provide  the  equi-energy 
density  contours  of  Figures  16  to  18.  These  contours  are  calculated  by 
joining  together  those  parallelepipeds  that  contain  the  same,  or  nearest, 
energy  value.  Each  has  been  formed  from  the  data  gathered  from  simulating 
20,000  incident  electrons  at  the  same  three  beam  voltages  as  those  of 
Figures  3  to  8  (10,  20  and  50  keV) ,  into  three  different  thicknesses  of 
PMMA  resist,  0.2,  0.5  and  1.0  microns.  For  example.  Figure  16  shows  four 
equi-energy  contours  for  a  beam  of  10  keV  electrons  when  incident  upon 
1.0,  0.5  and  0.2  microns  of  PMMA  respectively.  These  contours  are  drawn 
at  20%,  40%,  60%  and  80%  of  the  maximum  value.  Figures  17  and  18  show  the 
same  for  20  and  50  keV  beam  energies.  The  contour  energy  values  are  chosen 
quite  arbitrarily  and  it  is  intended  that  the  contours  should  show  the 
overall  energy  deposition  profile  and  no  attempt  is  made  here  to  relate 
these  to  a  developed  feature  profile  although  these  contours  must  be 
important  factors  in  the  shaping  of  the  final  feature  profiles.  The 
effect  of  using  a  thinner  resist  level  is  more  evident  from  these  diagrams, 
there  is  a  noticeable  improvement  in  resolution,  ie  less  broadening  of  the 
profile  towards  the  interface  of  the  resist  and  substrate,  in  the  0.2  micron 
resist  layer  over  the  0.5  and  1.0  micron  layers  in  both  the  10  keV  and 
20  keV  case  (Figures  16  and  17)  but  is  much  less  marked  in  the  higher 
50  keV  case  (Figure  18). 

Figure  19  shows  a  comparison  of  equi-energy  levels  within  a  0.5  micron 
layer  of  PMMA  after  being  exposed  by  20,000  electrons  at  the  three  different 
energies,  10,  20  and  50  keV.  Again  the  contour  energy  values  are  arbitrary 
and  consistent  through  the  three  diagrams  yet  it  is  clear  that  the  diagrams 
do  not  contain  the  same  number  of  contours.  The  10  keV  diagram  has  5  con¬ 
tours,  the  20  keV  diagram  has  4  having  lost  the  centre  contour  and  the 
50  keV  diagram  has  3  having  lost  the  inner  most  two.  This  is  because  the 
rate  of  loss  of  energy  by  an  electron  is  inversely  proportional  to  its 
energy  (see  equation  (4))  and  therefore  less  energy  has  been  deposited  by 
the  more  energetic  electrons,  hence,  the  resist  is  less  sensitive  to 
exposure  by  the  higher  energy  electrons. 


10 


Just  for  comparison,  the  equi-energy  contours  for  the  tri-level  resist 
structure  of  Figure  10  and  the  Gallium-Arsenide  case  of  Figure  9  are  pre¬ 
sented  in  Figures  20  and  21  respectively.  Comparing  these  contours  with 
those  for  the  0.5  pm  of  PMMA  on  silicon  structure  in  Figure  17  the  differ¬ 
ences  remarked  upon  earlier  become  clear.  Looking  at  Figure  20  the  resolu¬ 
tion  has  obviously  improved  somewhat  over  the  single-level  structure  and 
the  problems  of  increased  backscatter  through  using  a  more  dense  substrate 
are  very  well  portrayed  in  Figure  21. 

Figures  22-24  are  another  method  of  displaying  the  same  data  of 
Figures  16-18  and  show  the  normalised  energy  dissipation  distributions  at 
the  resist-substrate  interface  for  the  three  resist  thicknesses,  0.2,  0.5 
and  1.0  microns.  Figure  22  has  been  simulated  with  20,000  50  keV  electrons. 
Figure  23  with  20,000  20  keV  electrons  and  Figure  24  with  20,000  10  keV 
electrons.  The  problem  caused  by  backscattered  electrons  is  quite  evident 
in  these  diagrams.  One  notices  that  for  the  20  keV  and  10  keV  case 
(Figures  23  and  24)  there  is  a  non-trivial  amount  of  energy  dissipated  at 
one  micron  distant  from  the  incident  position,  especially  in  the  thicker 
resist.  This  will,  in  turn,  cause  inter-proximity  effect  between  neighbour¬ 
ing  features  that  are  less  than  1.0  microns  apart.  Further  to  these 
Figure  25  shows  the  dissipated  energy  distribution  in  1.0  micron  of  FMMA  on 
Silicon  at  five  equi-spaced  depths  into  the  resist.  Again  we  can  see  that 
even  at  the  surface  of  the  resist  there  is  a  background  exposure  that 
stretches  out  beyond  1.0  micron  from  the  point  of  incidence  and  this  will 
limit  the  maximum  exposure  contrast  achievable. 

The  dissipated  energy  data  can  also  be  presented  as  three  dimensional 
surfaces,  as  in  Figures  26  to  30.  Figures  26,  27  and  28  display  the  data 
of  Figures  17  and  23  showing  the  energy  deposited  by  a  20  keV  electron  beam 
of  spot  diameter  of  0.2  microns  in  1.0,  0.5  and  0.2  microns  of  PMMA  on 
silicon  respectively.  Figures  29  and  30  show  the  energy  deposited  by  a 
50  keV  and  a  10  keV  electron  beam  respectively,  of  spot  diameter  of  0.2 
microns  in  0.5  microns  of  PMMA  on  Silicon.  These  five  figures  are  included 
to  show  the  flexibility  of  the  data  achieved  from  the  Monte  Carlo  simulations 

No  energy  data  for  the  ion  beam  case  is  presented  here  because  it  is 
felt  that  secondary  ion  effects  need  to  be  taken  account  of  before  any 
serious  exposure  modelling  can  be  attempted.  The  primaries  trajectory 
plots  were  included  to  show  the  potential  for  high  resolution  lithography 
that  focussed  ion  beams  have. 

4.3  LIMITATIONS  OF  THE  MODEL 

No  account  is  made  for  the  mechanism  involved  in  the  exposure  of  the 
polymer  resist.  Adesida  et  al(?)  used  Gryzinski's  excitation  function  as  a 
replacement  for  Bethe’s  continuous  slowing  down  equation  to  try  and  model 
the  exposing  process  more  directly.  They  concluded  that  this  approach  gave 
a  more  realistic  answer  for  the  transmitted  electron  energy  distribution. 

The  effect  of  secondary  electrons  has  been  studied  by  Murata  et  al^8;  and 
again  they  produced  some  slightly  different  results.  The  main  effect  was  a 
broadening  of  the  energy  distribution,  especially  towards  the  surface  of  the 
resist.  However,  it  was  considered  that  the  difference  these  adjustments 
made  to  the  overall  results  did  not  warrant  the  increased  complexity  in  the 
models  as  this  is  primarily  a  study  of  the  lithographic  process  and  is  not 
an  in-depth  study  of  particle  scattering  in  solids.  A  further  limitation. 


within  the  Monte  Carlo  method,  is  that  both  substrate  and  resist  are  assumed 
to  be  totally  amorphous.  This  will  probably  have  a  bigger  effect  on  the 
electron  case  than  the  ion  because  any  channelling  due  to  the  crystalline 
structure  will  effect  the  backscatter  and  we  expect  zero,  or  very  nearly, 
backscatter  of  the  ions.  However,  we  are  only  interested  in  the  energy 
deposited  into  the  amorphous  resist  layer  in  this  application,  thus 
channelling  is  of  little  direct  relevance. 

The  Monte  Carlo  method  has  a  few  practical  disadvantages  such  as  the 
need  for  a  large  computer  and  hours  of  CPU  time  to  simulate  a  large  number 
of  electrons  in  order  to  smooth  out  the  statistical  fluctuations.  Also  the 
data  is  returned  in  histogram  form  instead  of  analytical  form.  However, 
these  problems  can  be  tolerated  because  the  model  is  fairly  simple  and  each 
run  of  the  programme  returns  enough  information  that  it  need  only  be  done 
once  for  each  substrate  and  resist  parameters.  This  data  can  then  be 
manipulated  further  for  changes  in  other  conditions,  like  beam  spot  size, 
afterwards. 

5  CONCLUSIONS 

The  Monte  Carlo  method  of  simulating  scattering  and  energy  loss  by  energetic 
particles  incident  upon  resist-substrate  structures  has  been  investigated  as  a 
means  to  understanding  the  exposure  of  resist  polymers  in  particle  beam  litho¬ 
graphy  . 

The  results  presented  clearly  show  the  problem  of  backscattered  electrons 
limiting  the  resolution  capabilities  of  electron  beam  lithography  and  suggest 
that  improvements  can  be  made  by  reducing  the  thickness  of  the  resist  layer 
employed  and  by  using  higher  energy  electrons.  The  trajectory  plots  of  the 
primary  ions  show  that  ion  beam  lithography  is  limited  only  by  the  small  amount 
of  forward  scattering  experienced  by  the  ions  and,  therefore,  has  the  potential 
for  high  resolution. 

The  models  used  in  this  study  are  comparatively  simple  yet  very  powerful  and 
are  able  to  indicate  the  roots  of  many  problems,  such  as  proximity  effect, 
experienced  by  electron  beam  lithographers.  The  models  can  be  used  to  investigate 
viable  solutions  to  these  problems  and  to  predict  performance  on  new  substrates 
and  new  resist  structures. 

6  FUTURE  WORK 

The  work  presented  here  is  currently  being  expanded  to  include  simulations 
of  other  processes  within  the  lithographic  system  to  improve  the  predictive  power 
of  the  simulations.  The  extensions  to  be  included  are: 

i.  Simulation  of  shaped  beam  systems.  This  is  a  very  simple  modification 
to  the  present  system.  All  that  is  required  is  that  the  data  produced 
by  the  Monte  Carlo  simulation  be  convoluted  with  the  appropriate  beam 
shape  required  (see  section  3.3). 

ii.  Inclusion  of  secondary  ion  effects.  As  mentioned  in  section  4.2 

secondary  ions  are  considered  to  play  an  important  role  in  exposure  by 
ion  beams.  Therefore,  these  will  be  included  in  the  simulations  for 
ion  beam  exposures. 
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iii.  Development  modelling.  To  complete  the  modelling  suite  multi-pass 

exposures  of  features  in  resist  and  their  realisation  by  development 
are  to  be  simulated.  This  will  provide  the  capability  of  being  able 
to  produce  theoretical  features  to  compare  with  those  we  gain  from 
our  experimental  programme. 
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Figure  9.  Simulated  trajectories  of  100  electrons  at 
20  keV  incident  upon  0.5  um  of  PMMA  on 
Gallium  Arsenide 


Figure  10.  Simulated  trajectories  of  100  electrons  at 

20  keV  incident  upon  a  tri-level  resist  'V-v-1 
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Figure  11.  Simulated  trajectories  of  50  lOOkeV  Gallium 
ions  incident  upon  0.5  microns  of  PMMA  on  a 
silicon  substrate. 


Figure  13.  Siaulated  trajectories  of  50  200keV  silicon 
Ions  incident  upon  0.2  aicrons  of  PMMA  on  a 
silicon  substrate. 


Figure  18.  Equi-energy  contours  simulat 
electrons.  The  energy  levels 
consistent  through  the  three 
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Figure  20.  Equi-energy  contours  within  the  top  resist  layer 
of  the  tri  level  resist  structure  of 
figure  10.  Simulated  with  20000  20  kev  electrons 
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Figure  21.  Equi-energy  contours  simulated  with  20000  20  keV 
electrons  into  0.5  um  of  PMMA  resist  on  a 
Gallium-Arsenide  substrate. 
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Figure  26.  Energy  deposition  surface  simulated  with  a 

0.2  um  beam  of  20  keV  electrons  incident  upon 
1.0  um  of  PUMA  on  Silicon. 
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Figure  27.  Energy  deposition  surface  simulated  with  a 

0.2  um  beam  of  20  keV  electrons  incident  upon 
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Figure  28.  Energy  deposition  surface  simulated  with  a 

0.2  urn  beam  of  20  keV  electrons  incident  upon 
0.2  urn  of  PMMA  on  Silicon. 
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Figure  29.  Energy  deposition  surface  simulated  with  a 
. . . „P»?  uw  beam  of  50  keV  electrons  incident  upon 
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